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Download HMI Magnetogram 
Line-of-Sight (LOS) or Vector 


MAG4 is 
Identify Active Regions LOS com pletely 
Use SHARPS for vector 
automated, from 
downloading 
magnetograms to 


= outputting and 


Predicts Event Rate storing forecast 


v = products. 


Makes Forecasts 


Calculates Free-Energy Proxy 








Ges, MAG4 Science Background S 





+ Flares and CMEs are known to be drivers of the most severe space 
weather 


+ Flares and CMEs typically originate in active regions (aka sunspots) 


+ Flares and CMEs are examples of exceptionally large explosive 
releases of magnetic energy stored in the corona 


+ While the amount of free energy cannot be measured directly, free- 
energy proxies can be measured 


+ Event rates have been shown to be correlated with the magnitude of 
the free-energy proxies 


ioe What is MAG4? S 





+ MAG4 (Magnetogram Forecast), developed originally for NASA/SRAG 
(Space Radiation Analysis Group), is an automated program that 
analyzes magnetograms from the HMI (Helioseismic and Magnetic 
Imager) instrument on NASA SDO (Solar Dynamics Observatory), and 
automatically predicts the event rate (or probability) of major flares (M- 
and X-class), Coronal Mass Ejections (CMEs), and Solar Energetic 
Particle Events 


+ MAG4 does not forecast that a flare will occur at a particular time in the 
next 24 or 48 hours; rather the probability of one occurring! 


* GONG (Global Oscillations Network Group) magnetograms, can be 
used instead as a backup but at a lower forecast accuracy 


+ Present cadence of new forecasts: 96 minutes (60 minutes for 
CCMC) Vector magnetogram actual cadence: 12 minutes 


42%% R20 Timeline of MAG4 S 





1973 The MSFC (Marshall Space Flight Center) vector magnetograph was built to support 
Skylab 

2000-present MSFC used vector magnetograms to study CME correlation with free- 
energy proxy 

2007-12 A DOD/Multidisciplinary University Research Initiative/Neutral Atmosphere 
Density Interdisciplinary Research helped support funding of the basic research 

2008 Partnered with JSC/SRAG (Space Radiation Analysis Group) and won an R20 
NASA/Technical Excellence Initiative grant: Began building a database that grew to 
~40,000 magnetograms of ~1,300 active region, covering years 1996-2004 with event 
catalog from SOHO/MDI (Solar and Heliospheric Observatory/Michelson Doppler Imager) 
observations 

2010-present NASA's HEOMD (Human Exploration and Operations Mission Directorate) 
support 

2010 SDO is launched began transitioning from MDI to HMI line-of-sight magnetograms. 
2011 MAG4 installed at SRAG a NRT (Near-Real-Time) forecasting tool, and SRAG 
began pre-operations testing 

2012 Provided NOAA web access to MAG4 NAT forecasts 

2013 Improve MAG4 so that it can use a combination of free-energy proxy and previous 
flare activity 

2015 Transition to HMI line-of-sight to vector magnetograms 

~AMS2016 





Calculating the Free-Energy Proxy S 





« When the transverse A magnetogram of an active region 
gradient of the vertical (or 
line-of-sight) magnetic field is 
large, there is more free- 
energy stored in the 
magnetic field 


+ For each Active Region: 
The integral of the gradient 
along the neutral line is the 
free-energy proxy 
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Converting Free-Energy Proxy to Predicted Event Rates S 
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Forecast Curves Free-Energy Proxy W(kG) 


These empirical forecast curves are used to convert our free-energy proxy 
into predicted event rates. Curves are derived from a sample of 40,000 
magnetograms, from 1300 active regions observed between 1996-2004. 


~AMS20i6 OSS 


+ Magnetograms are spatial 
maps of the magnetic field ss pene e eye 
strengths 7 


* They come in two basic types 
line-of-sight (right) 
- vector magnetograms 


+ Free-energy proxies can be 
measured for Active Regions 
(areas with sunspots) from 
either type of magnetogram 


+ Line-of-sight magnetograms 
suffer reduced accuracy 
further from disk center 


A full-disk line-of-sight magnetogram of 
roel) 220) 


Magnetograms & identify ARs 
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Forecast 


Multiplicative uncertainty example 


Rate 1 sigma Probability 
Events/day 66% Confidence 

0.02 0.01-0.05 0.7-5% 
0.7 0.3-1.9 20-80% 


For a Multiplicative Uncertainty of 2.7x 


MAGA Forecast 








2012/03/06 14:23 
# ARH WL, Lng Lat 


(Ke) (dea) 
3148 8 et -7 
511420 650-1 17 
611490 11-25 20 
431186 -27 


Disk Forecast Rates 
Multiplicative Uncertainties 
Disk All-Clear Forecast Probabilities 
Uncertainties 





24 Hour Event Rate Dist 
M&X CME FCME X SPE (deg) 
0.020 0.020 0,009 0.002 0.003 27 
0.700 0.400 0.200 0.100 0.080 44! 
0.020 0.030 0.010 0.004 0.005 32! 
0.000 0.001 0.001 0.000 0.000 45! 
0.800 0.400 0.200 0.100 0.090 

27x 21x 22x 30x 24x 

50.00% 70.00% 80.00% 90.00% 92.00% 
40.00% 20.00% 10.00% 10.00% 7.00% 


f Safe and Not Safe Days 


June 26, 2013 March 7, 2012 
C1, C1.5 flares X5.4, X1.3, C1.6 
CME 2684, 1825 km/sec, 


Solar Energetic Proton Event reaches 
6530 particle flux unit >10MeV 











‘ley How Well Does MAGA Forecast: S 








2. Skill Metrics Significance of Upgrade 


Mcintosh/NOAA 259 638 631 18476 93.7 0.29 0.71 0.26 0.26 


Free-Energy Proxy 273 284 618 18830 95.5 0.31 0.50 0.35 0.47 
Present MAG4 


Free-energy proxyand 340 317 551 18797 95.7 0.38 0.48 0.42 0.49 
previous flare activity 
Upgraded MAG4 


Best 890 0 0 19114 100 1 0 1 1 


Improvement in Metric 


Mcintosh/NOAA 1.840.5 0.0340.05 0.2140.07 0.10+0.04 0.21+0.07 
Present MAG4 (40) (0.30) (30) (20) (30) 
Present MAG4 0.240.2 0.0840.03 0.02+0.05 0.06+0.03 0.03+0.05 
Upgraded MAG4 (0.70) (20) (0.50) (20) (0.50) 






Rote of No X or M flore 





How Well Does MAGS Forecast: S 





During periods when flare-productive active regions cross the disk, the 
predicted rate and actual rate both increase, providing situational awareness 
The results are best when flares and predicted rates are limited to inner 45 
degree circle (Right) 

M or X-Class Flares 


Full disk Rates Central Disk 45° 
5 Day smoothed X&M flare rate NG Day smoothed X&M flare rate 


Predicted Event Rate 


Rote of No X or M flare 
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Gx, MAG4 Improvements: Vector Magnetograms 









* Both vectors shown in 
red have positive B, 
(magnetic field out of Actual Examples 















the sun), but have 
opposite sign B, gg and False Neutral 
thus a false Lines occur on 


(unphysical) neutral line limbward sides of 










in the line-of-sight Se DO 

(LOS) field. 
Problem fixed by 
converting from 
Bios and 


Briansverse 10 Bz 
and Byorizontal 






+ For LOS most accurate to 30 
heliocentric degrees. (RED) : Peer eye) 


+ For Vector similar accuracies 
should exist to 60 heliocentric 
degrees. 


A full-disk line-of-sight magnetogram of 
rool) 220) 


Magnetograms & identify ARs 





‘AMS 2076 





Near Future Improvements S 





+ MAG4 presently AGS by NASA-MSFC/ 

— Predicts what MDI 
will measure if AR 
was at disk center. 

— Then uses the 
forecast curves from 
MDI to forecast 
events. 

+ MAG4 needs to 
upgrade to take full 
advantage of vector 
magnetograms. 


*  htto:/Awww.uah.edu/csp 
ar/research/mag4-page 
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is) How Well Does MAG4 Forecast: S 


J ctual ves ActualNo | 


Predict Yes YY YN 
Predict No NY NN 


Metric Equations 


Percent Correct PC=(YY+NN)/(YY+YN+NY+YY) 

Probability of Detection © POD=YY/(YY+NY) 

False Alarm Rate FAR=YN/(YY+YN) 

Heidke Skill Score HSS=2*(YY*NN-YN*NY)/[(YY+NY)* 
(NY+NN)+(YY+YN)*(YN+NN)] 

True Skill Score TSS=(YY*NN-NY*YN)/((YY+NY)*(YN+NN)) 


Ce. How Well Does MAG4 Forecast: S 





PC Percent Correct (0 to 100%) 
ital eta POD Probability of Detection (0 to 1) 
Predict Yes YY YN FAR False Alarm Rate (1 to 0) 
Predict No NY NN HSS __Heidke Skill Score (-1 to 1) 


TSS True Skill Score (-1 to 1) 


Frorcastnethod YY [YN [NY [NN [PC(%) [POD] FAR [HSS [TSS | 


MciIntosh/NOAA 259 638 631 18476 93.7 0.29 0.71 0.26 0.26 
Free-Energy Proxy 273 «284 «©6618 §=18830 95.5 0.31 0.50 0.35 0.47 


Free-energy proxyand 340 317 551 18797 95.7 0.38 0.48 0.42 0.49 
previous flare activity 


Best 890 0 0 19114 100 1 0 ul 1 
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Ge, Improvements with Vector MAG4 S 





+ Correct for systemic errors in measurements due to radial dependence 
of transverse field noise and foreshortening. 


* Access skill scores relative to present LOS 

+ Investigate if a vector measure is better than LOS measure 
+ Investigate if change of measure is important 

+ Investigate if combination produces better measure 

* Correct SPE forecasts based on Longitude of source region 


+ Partner with heliosphere modeling to better forecast where the particles 
would go. 


* Twin CME 





Ges, Basic Research to do for later R20 








Determine if growth/decay phase 
of an active region has additional 
prediction ability 


Determine if looking at twin-CME 
or single-CME can produce 
improved SEP forecasts 


“Wlse (G) 

















Ma 
4 5 6 7 8 99 10 11 12 13 14 
November 2004 November 2001 








Flore Mognitude 









Research to Operations and Operations for FY14 S 


& 


+ Continue to supply support, and upgrades to NASA/SRAG 





+ Finish and evaluate using vector magnetograms 


+ Fully implement in consultation with SRAG using combined free-energy 
proxy and previous flare activity 





+ Improve All Clear Forecast for SEP 1000.0 


+ Partner with Air Force, and through them NOAA 100.0 
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Improving the Forecasts: S 





X&M Flore Rate (day~') 


1. Recent Flare History 














1.000 = 1,000 
> 
3 
0.100 3 0.100 
3 — Free Energy Only 
0.010 z 0.010 —— Recently Flaring 
a — Recently Non-flaring 
0,001 0,001 
1 10 100 1 10 100 


Free-Energy Proxy W(kG) Free-Energy Proxy W(kG) 


Active regions that have recently produced an X- or M-Class flare are 
more likely to produce flares in the near future 





Ge. Free-Magnetic Energy S 





Contours Vertical Magnetic Field 


Arrows Transverse Magnetic Field Currents ~10'? Amps 
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slxy How Free Energy Proxy Evolves and Clustering of Flares S 
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Free-Energy Proxy W(kG) 
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Evolution of Free-Energy Proxy 
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24—Hour X&M Flare Rate 


8 9 10 #11 
March 2012 


Forecasted X&M Flare Rate 









X-Class Flares * 
M-Class Flares + 





6 7 8 9g TOs: 01 
March 2012 


The free-energy proxy evolves on time periods of days => 
the forecast is on those time scales 





AMS 20iS 
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Flares occur when high free-energy proxy active regions crossing disk. 


How well it works (All-Clear) S 





All Clear Probability 


o 
2 

2 

= 

& 

° 

x< 

2° 

Zz 

‘o 

2 

2 

2 Times of X or 
& M Flares 





O9—Apr 02-—May 26—May 18—Jun 
Date 2012 





